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ABSTRACT 

The origin of complex organic molecules such as amino acids and their precur- 
sors found in meteorites and comets is unknown. Previous studies have accounted 
for the complex organic inventory of the Solar System by aqueous chemistry on 
warm meteoritic parent bodies, or by accretion of organics formed in the inter- 
stellar medium. This paper proposes a third possibility: that complex organics 
were created in situ by ultraviolet light from nearby O/B stars irradiating ices 
already in the Sun's protoplanetary disk. If the Sun was born in a dense clus- 
ter near UV-bright stars, the flux hitting the disk from external stars could be 
many orders of magnitude higher than that from the Sun alone. Such photolysis 
of ices in the laboratory can rapidly produce amino acid precursors and other 
complex organic molecules. I present a simple model coupling grain growth and 
UV exposure in a young circumstellar disk. It is shown that the production may 
be sufficient to create the Solar System's entire complex organic inventory within 
10 6 yr. Subsequent aqueous alteration on meteoritic parent bodies is not ruled 
out. 



Introduction 



The early Solar System was replete with complex organic molecules, which can be seen 
today in preserved ancient bodies such as meteorites and comets. Upwards of 100 different 
amino acids have been detected on chondritic meteorites such as Mur chison, Allende anc 



Orguie l. Many of these have no known natural terrestrial occurrence (jEhrenfreund et al 
200 lbl ). and are believed to be of extra-terrestrial origin. Comets and the interstellar 
medium (IS M) are also r ich inventories of complex organic molecules, including amino acid 



precursors (I Snyder 



20061 ). Amino acids and other complex organics have have implications 
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for the origins of life on Earth, so understanding their formation and history is an area of 
great interest. If these molecules are formed and distributed easily in a variety of disks 
and conditions, then pre-biotic compounds may be common throughout distant planetary 
systems. 

Of the complex organics in the solar system, amino acids are particularly interesting 
to study. These compounds are necessary for life, and may have been biotic precursors 
on Earth. Carbonaceous chondrites are typically 5% or more carbon by mass, most 
o f which is in aromatic poly mers and thousands of other pre-biotic organic compounds 



(T IHayatsu and Anders 



198ll ). Laboratory stud ies have found Murchiso n to have some 



19901 ) . The simplest 



10-30 ppm by mass in identified amino acids (jShock and Schultd 
amino acid (glycine, NH^CH^ COOH) has recently been detected in comets, bu t so far has 



eluded detection in the ISM ( lElsila et al. 



2009; 



Snyder 



2006; 



Kuan et al 



20Q3J ). Potential 



amino precursors su ch as formic acid ( HCOOH) have been detected in comets but not the 



interstellar medium (jHollis et al. 



2003 



Bockelee-Morvan et al. 



2000J). However, amino acids 



in ice form in the interstellar medium (ISM) would be more stable but harder to detect 



than their corresponding gas forms in comets (jEhrenfreund et al. 



2001a 



Chyba et al. 



1990), 



so it is likely that they exist in the ISM but have not yet been detected. 

The origin of these molecules is unknown. Two broad explanations exist for the 
amino acids in our solar system. First, they may be produced endogenically, by chemical 
synthesis within the solar system itself. A variety of energy sources for this exist, including 
infall heating, radiogenic heating, lightning, and shocks. For the amino acids present on 
meteoritic parent bodies, the best-studied endogenic process is Strecker synthesis. This is 
a method by which amino acids are formed in aqueous environments, such as the warm 



sub-surface aquifers that could have been present o n asteroids as they were 



and other radionuclides during their first few Myr ( jEhrenfreund et al 



reated by 26 A1 



2001bJ). Amino acids 
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are produced in the laboratory by th is process, and suffici ent 26 Al existed to heat parent 



bodies to liquid water temperatures ( IMcSween et al. 



20021 ). However, a serious problem 



exists with this model. Isotopic measurements of the chondritic amino acids consistently 
show deuterium enhancements of SD = 600-2000 %o, while measurements of the water in 
these same bodies show deuterium depletions of roughly 100 %o- In situ synthesis of aminos 



does not preferentially chan ge the 



come from distinct sources (ILerner 



2 /H ratio, so this argues that the water and organics 



19971 ) 



Second, the molecules may have been inherited exogenically from the interstellar 
medium which formed the solar nebula. The ISM is known to be rich in complex 
chemistry, with over 150 differe nt gas-phase species detected to date in molecular clouds 



(IHerbst and van Dishoeck 



20091 ). Amino acids have not been detected in the ISM, but 
sugars, alcohols, polycyclic aromatic hydrocarbons (PAHs), and other complex molecules 
have been. The formation of these compounds in the ISM is thought to be due to a 
combination of processes, including gas phase, gas-grain, and UV-grain reactions. Once 
molecules are formed, stable compounds can be incorporated into young disks as new 
YSOs condense out of the ISM. This pathway is supported by m easurements showing th e 



organic composition of comets and the ISM to be quite similar (IBockelee-Morvan et al. 



2000 



Irvine et al. 



2000) 



Laboratory results have shown that simple ices can be turned into a zoo of amino 
acids and other complex organ i cs simply by the presence of UV flux followed b y a warming 



stage (INuevo et al. 



2008 



2007 



Bernstein et al. 



2002 



Munoz Caro et al 



20021 ). Though 



these studies are performed at fluxes higher than the ISM, their total photon dosages 
are comparable to the current problem. Amino acids similar to those seen in meteorites 
are created. Moreover, chemical reactions (including photolysis) which occur at cold 
temperatures < 70 K can preferentially increase D/H because deuterium's higher mass 
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allows it to bond more readily than H at low temperatures. D/H enrichments are seen 
throughout the ISM, and the D/H enrichments in amino acids suggests that they too have a 



low-t emperature origin compared with the water in meteorite parent bodies (jSandford et al 



200ll ) . Continued UV exposure may sometimes destroy the same large molecules the UV 
created earlier, but for a proper range of timescales and fluxes, cold-temperature photolysis 
of ices may be a more plausible pathway to forming the meteoritic amino acids than warm 
aqueous synthesis. 

The ISM formation model is not without problems. First, the high optical depth 
within dense molecular clouds blocks nearly all external UV light. The only source of 
UV within the clouds is that caused by cosmic rays interacting with gas in the clouds, 
resulting in a UV flux of 10 ~ 4 -1CU 5 Go, where Gq is the interstellar flux at the Sun today 



( jPrasad and Tarafdar 



19831 ) . Second, some organic molecules, such as the CHON particles 



in comets like Halley, may be easily destroyed by shocks during infall, although some 



deba te exists on this point fjVisser et al. 



2009 



Mumma et al 



1993 



Zahnle and Grinspoon 



19901 ). Finally, some of the common meteoritic amino acids, such as a-aminoisobutyric 



acid (AIB) and isov aline, have not yet been produced by irradiation in the laboratory 



([Hudson et al. 



20081 ). and it is not known whether th is is due to these species 



Ehrenfreund et al. 



only being 



2001bl ). or simply 



formed during Strecker synthesis (as hypothesized by 
the ice experiments not using the proper initial mixtures. Nevertheless, overall the ISM 
formation model paints a plausible picture and may well play a role in the solar system's 
organic history. 

This paper proposes a third pathway for the formation of these complex organic 
molecules, which has not been examined previously. In our model, the solar nebula forms 
within a large molecular cloud similar to those in Orion. The Sun and its disk form 
completely, and condense and begin to grow. Within the next 1 Myr, nearby O and B 
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stars turn on, bathing the young disk in UV light. These UV photons photo-evaporate gas 



from the disk ( IThroop and Bally 



2005 



Johnstone et al 



19981 ) . but also irradiate small ice 



grains exposed in the disk's outer skin layer. The simple ices grains, such as H 2 0, CO2 
and NH3, are exposed to the UV and begin to photolyze into more complex species. Each 
individual grain is exposed only briefly to UV light; they spend most of their time in the 
disk's dark, turbulent interior. Grains continue to grow, gas which has not formed planets 
is lost due to photo-evaporation and viscous loss, and within 5 Myr a gas-free debris disk 
is left with its ices enriched in complex organics. These organics can then be incorporated 
into comets, asteroids, and the planets. Organics produced in this way could complement, 
and perhaps greatly exceed, those produced by the other two pathways. This method is 
similar to the exogenic ISM production in that it relies on UV photochemistry, but at far 
high flux (10 6 Go vs. 10~ 5 Go), at warmer temperatures, for a much shorter time. The 
model allows for subsequent aqueous alteration as seen in the meteoritic record. 



Work by 



Robert 



(J2002J) proposed that the organic material was created by X-ray 
irradiation of the disk by the Sun during its T Tauri phase. They did not present a disk 
model, but used high-precision D/H measurements of organic and non-organic material 
to show that D/H fractionation varied with heliocentric distance, as would irradiation. 



Remusat et al. 



(120061 ) compared the meteoritic and interstellar D/H values, along with their 
C-H bond dissociation energies, and concluded that the solar system's D/H enrichment was 
created in situ, rather than inherited from the ISM. They proposed that the young Sun 
might provide the necessary UV source; their model did not study the timescales or fluxes 
involved. 

The source of the Earth's organic inventory (as opposed to the Solar nebula's) 
is a parallel question which has received some attention. It is believed that although 
some organics were probably synthesized on the young Earth in situ by lightning or 
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spark discharges ([M iller 



( IChyba and Sagan 



and Ureyl Il959f ). sho ck heating in the terrestrial atmosphere 



19921 ) . warm ocean vents (1 Corliss 



1990), or any number of other 



terrestrial processes, far more were probably delivered from external s ources such as comets 



asteroids, and interplanetary dust particles ( IChyba and Sagan 



19921 ). This paper does 



not further address the origin of life or the Earth's inventory, except to acknowledge that 
increasing global abundance of organics in the solar nebula probably results in increased 
delivery to the Earth as well. 

This paper uses a simple model to describe the production of complex organics in 
UV-illuminated disks in a variety of cases. The results are necessarily general, and do not 
explain the abundances or species in one particular sample or disk, but provides an initial 
assessment of the problem. External UV photolysis has been ignored in almost all previous 
models, yet it may be one of the most important sources of energy in both the young 
Solar System and other proto-planetary disks. The problem is set up in SJ21 and the model 
described in §31 §H contains results, which are discussed in §5j Conclusions are in §61 and 
Appendix A contains a derivation of the simple grain growth model used here. 



2. Background 



2.1. Energy sources in the early solar nebula 



Various surveys have shown that the majority of st ars are born in dense clusters of 
300- 10 4 stars, where massive O and/or B stars can form (lAdams et al. 



2006 



Lada and Lada 



20031 ). Due to the presence of decay products from short-live radionuclides such as 60 Fe that 
are only produced in supernovae, our Sun is th ought to have formed in such an environment 



(jHester et al. 



2004 



Tachibana and Huss 



20031 1 . The effects of such an environment can 



be enormous, compared with the classical 'closed box' view of Solar System formation 
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where any environmental effects are easily ignor ed. Photo-evaporation by external stars can 



remove the disks on Myr tim escales or shorter (IThroop and Bally 



2003 



Johnstone et al. 



in the outer disk ( iDuncan et al. 



2005 



Matsuyama et al. 



1998 LC 



ose approaches between s tars may be important, especially 



2008 



Adams et al. 



20061 ). The high gas densities in such 



clusters can cause late infall of mo 



begun to form (IThroop and Bally 



ecula r cloud material onto disks after planetesimals have 



20081 ). And, accretion of ejecta from supernovae and 



massi ye stellar winds may contaminate the disk, changing its structure, composition, or 



both ( IThroop and Bally 



2010 



Cunha et al. 



19981 ) 



The effect of the local environment on solid-phase chemistry in t he Sun's pre-p l aneta ry 
disk has rec e ived s c ant attention. P assing r eferen c es can be fou nd in 



Irvine et al. 



(11993 ) . 



(2000), 



Throopl (120001) . iFegleyl (119991 ). iPrinnl (11993 ). and 



Bergin et al. 


( 


2007) 


van Dishoeck et al. 



Prinn and Fegleyi ( 119891 ) provides a complete review but is now somewhat eclipsed 
by more newer RT models in the astrophysical literature (discussed below). Their paper 
does not directly address the formation of organics, but examines the energy budget for 
chemistry from various sources in the young solar nebula. They find that the largest 
energy source is the gravitational energy given off by the nebula as it collapses. In addition 
to the entire bulk disk he ating, individual grains are heated an d may sublimate ices as 



they fall into the nebula ( iVisser et al. 



2009 



Lunine et al. 



19911 ) . While this may cause 



some fractionation and possible hydrolysis of organic molecules, simply warming ices from 
10 K to 200 K for a few hundred seconds during the infall is not itself a pathway to form 
amino acids. However, their second-largest source, shocks and lightning, certainly is. They 
estimate the total usable energy for reactions by these processes to be roughly 4 x 10~ 6 
that of the total gravitational collapse energy of the disk. Radionuclides {e.g., 26 Al) were 
estimated to provide a few orders of magnitude again less energy. They also estimated 
energy from photochemistry, based on fluxes from both the Sun and interstellar sources. 
They argued that photochemistry in the solar nebula was unimportant, especially from from 
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internal sources. However, more sophisticated recent models show that photochemistry can 
in fact be the dominant energy source. 

Internal UV sources. Several recent models for disk formation consider UV radiative 
transfer from the central star to the inner disk. The model of 



Woitke et al. 



(120091 ) handles 



photochemistry using a Monte Carlo method in a 2D disk of mass 1 MMSN, surrounding 
a T Tauri star. Over 70 chemical species are considered, including five ices, with a total 



of nearly 100 different chemical re actions. A broa dly similar model with 



described by lJonkheid et al. 



ewer species is 



(120071 ). The model by iGorti and Hollenbachl (120041 ) describes 
the UV photochemistry of 73 gas species in a 10 Myr old debris disk, putting it in a much 
lower optical depth regime than the other models. All three models are computed for a fixed 
age and do not evolve. The grain size distribution in all is uniform across the disk. All three 
models present sophisticated pictures of the radiative transfer and photochemistry within 
the disks, allowing predictions to be made for abundances and line strengths. However, the 
model does not include ice photolysis, grain growth, differences in grain size across the disk, 
or external illumination. The UV field reaches ~ 10 4 Go inward of 10 AU, but this is still 
much lower flux than the 10 6 Gq or more from external illumination. 



In calculating the energy budget of the nascent solar system, iPrinn and Fegleyl (119891 ) 



estimated that the the line-of-sight UV optical depth through the disk to be "7 million to 
110 million orders of magnitude." In their model, this extreme optical depth would prevent 



solar-driven photochemistry virtually anywhere outside the central 0.35 AU. _ E 
assumed a direct line-of-sight from the source was required, which 



oweyer, they 



Gladstone! ( 119931 ) points 



out is not the case because solar Lya photons can reflect off of interplanetary hydrogen far 
above the disk plane, thus illuminating the disk by reflected solar light. Considering the 
higher density and reflectivity of the interplanetary medium (1 PM) at the time, and the 



increase in early solar UV flux at 10 4 x or so over present values, 



Gladstone! (119931 ) calculate 
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that the young solar system was exposed to a solar UV flux roughly 10 times the present 



solar value. Scattering of solar Lya photons off of wi nds, jets and infa l 



Hollenbach et al 



ing m aterial can 



(11994 ) showed that 



also be a large source of UV onto the disk. Similarly, 
young stars can create a thin atmosphere above their disk, and this disk provides a source 
of indirect UV illumination which can drive photo-evaporation (and thu s photolysis), at 



rates 1 0-100 times lower than e xternal illumination would, or 10 4 -10 5 Go ( iMatsuyama et al. 

3 



2003|). 



Alexander et al 



fl2006al 



j) also looked at the case of photo-ev aporation from the 



centra l star in flared disks, which avoids the line-of-sig ht problem in 



Prinn and Feglev 



Alexander et al 



fl2006bl ) calculated 



( 119891 ). Combined with loss to inward viscous evolution, 
disk loss timescales of several Myr, consistent with observation of disk dispersal in young 
T Tauri stars. None of the papers by Gladstone, Hollenbach or Alexander considered the 
photochemistry effects of UV, but between these results and the detailed RT models above, 
it becomes clear that photochemistry from the central star can be important. Table [3] lists 
approximate fluxes for the various sources. 

External UV sources. While the interstellar UV flux in Taurus-like regions is low 
(~ 1 Go outside the cloud, and 10 -4 Go inside), stars in dense clusters such as the Orion 
Nebula Cluster (ONC) are subject to fluxes from O and B stars on the order ~ 10 6 Go- Ices 
in these clusters' disks are exposed to high levels of UV irradiation, enhancing the disks' 
abundance of organic molecules. In contrast to the internal sources, the external sources 
are stronger by 10 2 — 10 4 x, and illuminate the entire disk evenly rather than dropping off 
with radial distance. 

This paper examines only the effects of the external flux. Existing models provide 
sophisticated treatments of the internal UV flux and the gas photochemistry; our goal 
here is to understand the broad global effects of external flux on solid-state chemistry, in 
preparation for a more detailed model. 
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2.2. Energy budget 

A quick calculation shows the importance of external irradiation relative to other 
energy sources. First, consider the flux intercepted by a disk from its central star of 
luminosity L. In a disk with a radius:half- height ratio of 10:1, the fraction of flux intercepted 
by the disk is ~ 1/10, and the total energy deposited at all wavelengths in time At is just 

E Q ~ —L At. (1) 

° 10 V ; 

The total amount of energy available for chemistry through thermal heating, shocks, 
lightning, and so forth can be no greater than the disk's total collapse energy from the 
cloud to disk inner radius Rq, given a disk mass M d and stellar mass M s , or 



GM d M s 

E c = — jf- 1 - (2) 
Finally, the total UV energy absorbed by the disk from external stars is roughly 



£ uv = \T(2)7rR 2 d F ext At. 



(3) 



Typical values for solar-mass stars in dense clusters are are M s = 1M Q , M d = O.O1M , 
L = L Q , R d = 100 AU, R = 5 AU, F ext = 10 6 G , and At = 1 Myr. Plugging in, we find 
the ratio E so \ : E c : E\jy is approximately 1 : 100 : 10 4 . That is, in 1 Myr, the external UV 
dose received by the disk exceeds by 100 x the entire direct energy input from the central 
star, and the external UV dose exceeds by 10 4 x the entire collapse energy of the disk. Q 



Moreover, although therma 



reactions ( iPrinn and Fegley 



gradi ents can limit the amount of energy available for chemical 



19891 ). the UV energy is deposited directly into ice and dust 



l \{ the opposite were the case, photo-evaporation could not remove the disk. 
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grains, where individual photons cause photolysis. Therefore, the disk's dominant energy 
source is external flux, and this energy source's effect on chemistry cannot be ignored. 



2.3. UV Photolysis of Ices 



The experiments of iMillerl (119531 see also lUreyl 119521 ) showed that electric discharges 
within an atmosphere can rapidly create amino acids. Since this time similar experiments 
have been repeated with different initial species, temperatures, phases, and energy sources, 
finding the same general result that amino acids are relatively easy to produce given 
sufficient energy. In an astrophysical context, amino acids can be produced with energy 



(Bernstein et al. 


2002; 


Chvba and Saean 


1992: 


Miller anc 




3ada 


1988: 


Barak and Bar- Nun 


1975 




Miller and Urev 


19591). According to Hudson et al. 


( 


2008 


), "one conclusion is that 



energetic processing of almost any organic ice that contains C, H, N, and O probably 
results in the formation of amino acid precursors, which can be hydrolyzed to give the acids 
themselves." 

Most interesting in the regions of cold space are laboratory experiments in the last 
decade involvin g UV irradi ation of ice mixtur e s. Fo r instance, experiments by both 



Bernstein et al. 



(120021 ) and 



Munoz Caro et al. 



( 120021 ) started with thin ice mixtures 



including H 2 0, CH 3 OH, and NH 3 , chosen to simulate ISM conditions. The micron-thick 
ice mixtures were formed at 15 K, and then irradiated for 30 minutes at 10 7 Go- After 
warming, analysis in a gas chromatograph detected six amino acids, including glycine and 
alanine. Their laboratory irradiation corresponded to 10 7 yr in the interior of a dense 
cloud (10~ 4 G ), or 12 hours in the Orion nebula (10 6 G ). The net photolysis efficiency 
was around 0.25% (complex organics produced per photon), and by the end of their runs 



the majority of their C and N had been converted into new compounds (IBernstein et al. 
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19951 ). More recent experiments by lNuevo et all ( 120081 . 120071 ) measured upwards of a dozen 
different amino acids formed from irradiation of H 2 0, CO2, and NH 3 . 

The authors of all three experiments hypothesized that long exposures of cold ices to 
UV light could explain the multitude of organic molecules found in interstellar regions. If the 
solar nebula collapsed from such an enriched region of the ISM, subsequent incorporation 
of these species into preserved primitive bodies such as Murchison would be a natural 
consequence of their birth environment. 

Although UV can lead to the synthesis of organic molecules, excessive UV flux can 
be damag ing to these same mole cules, as evidenced by terrestrial sales of skin-protection 



products. 



Bernstein et al. 



( 12004 ) looked into the UV destruction of amino acids, and found 



that given enough flux they eventually formed nitriles, which are perhaps an order of 



magnitude more stable against further destruction. Similarly, 



Ehrenfreund et al 



( I2001ah 



found tha t several hundred yea rs of illu mination at 1 Go destroye d most of the aminos in ice 



matrices. 



Hudson et al. 



(12008 



see also 



Garrod and Herbst 



20061 ) showed that irradiation 



and hydrolysis of nitrile-containing ices results in amino acids being created again, so 
long-term irradiation may eventually result in an equilibrium between aminos and nitriles. 

In Orion, the flux hitting a disk varies with time as a star orbits through the cluster, 
changing its distance and tilt angle to the external O/B stars. Typical distances from 
the Orion core rang e from about 1 pc to 0.01 pc, causing the impinging flux to be in th e 



range 10 3 - 10 7 Go (IFatuzzo and Adams 



2008 



Throop and Bally 



2008 



Bally et al. 



19981). 



10 6 Go represents the flux hitting the disk skin. However, young disks are very optically 
thick, so the vast majority of the ice grains are nearer the midplane and thus sheltered from 



UV. Simple calculations by 



Throop 



( 120001 ) based on the disk mass and small grains gave 
typical initial optical thicknesses of close to 10 6 at 10 AU; thus, these individual grains go 
through a "broil-cool, broil-cool" cycle as they circulate through the disk, only occasionally 
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being exposed to the impinging flux. Disk midplane temperatures can be as warm as 50K 



out to 100 AU, so any nitriles produced during the i rradiation are like 



amino acids during w hen they shaded but still warm (IVisser et al. 



2009 



y to go on to 



Watson et al 



orm 



D'Alessio et al. 



2006|). 



20071 : 



Disk Model 



Throop and Bally 



he simple model used here is a modified version of that presented in 
(120051 ). The current model is ID and considers grain growth and external irradiation as 
a function of orbital distance R from a 1 M star. The disk is azimuthally and vertically 
homogeneous, with initial parameters specified in Table [TJ The initial disk mass is 
0.04 Mq. The disk is made of gas and dust (ice), in an initial mass ratio of 100:1. The 
disk has 40 logarithmically spaced radial bins, and it is evolved for 3 x 10 5 yr using a 
self-adjusting time-stepper coupled with a Crank-Nicolson integrator. At the end of this 
time, photo-evaporation has largely removed the gas disk. Real disks might continue to be 
exposed to UV for several Myr, but because of rapid grain growth, most of the irradiation 
of small grains happens at the beginning of the simulation and the results are not strongly 
sensitive to the time cutoff. 



3.1. UV Photolysis 



UV photons can create complex organic molecules or their precursors from simple ices. 
The effect can be roughly parametrized with a 'photolysis efficiency' e p , where p ~ 0.25% in 
molecules per photon for conversion from simple ices into amino acids. This parametrization 
skips the complex che mical chains involved, but is roughly the correct yield from pure ice 



(IBernstein et al. 



2002). 
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At each timestep, the model computes the UV exposure onto the disk. The 
instantaneous UV flux is 

F (R,t) = ^(l-exp(-r(R,t))) (4) 
where the optical thickness t(R, t) is 

assuming particles of size r and density p. The surface density S rf is that of ice grains alone, 
and does not change during the simulation. We assume scattering efficiency Q sca = 1, which 
is appropriate for all but the very smallest grains and does not change the results. 

At each timestep the flux F(R, t) is computed, and the total number of photons 
intercepted by each bin is recorded. These are then used to compute the total UV exposure 
in photons per molecule, 

* Wi ) = £^)J-A f , (6) 

where m m = 18 amu is the typical molecular mass. 

The present model assumes a fixed external flux of 10 6 G . This flux is what the 
well-studied Orion proplyds at 0.1 pc receive, and is picked for easy comparison with 
previous results on those disks. Two additional cases look at higher and lower fluxes. 



3.2. Grain Growth 



Grain growth is important to the model, because once surface area is locked up in 
large grains, it cannot be exposed to UV radiation. Both observational and theoretical 
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arguments shows th at grains grow rapid 



as short as decades (jDominik et al. 



2007 



y, reaching centimeter or meter sizes on ti mescales 



Johansen et al. 



2007 



Throop et al. 



200 lh . Grain 



growth (together with gas loss) ca uses the disks to clear in 



out, typically within 1-5 Myr (e.g. 



optica l depth from the inside 



Alexander and Armitage 



20071 . and references therein). 



The physical mechanisms for the early stages of grain growth are not well u nderstood. 



However, a simple parametrization is useful for the present study. 



Throop! ( 120001 ) performed 



semi-analytic calculations of accretionary grain growth of a size distribution n(r,R,t) of 



particles. For grains up to a few cm, the grai ns were coup 



velocities depending on the gas eddy speeds of 



Mizuno et al 



e d to the g as with colli sion 



(|1988h and 



Voelk et al. 



(1980) 



while for larger sizes they used particle-in-a-box collision probabilities. In both regimes they 
assumed a constant sticking probability e s . After integrating their distributions numerically, 
they combined an analytic expression for grain growth with numerically-determined 
coefficients from the eddy velocity simulations to describe a the net grain growth. In their 
model, adopted here and derived in more details in Appendix A, the typical grain size r 



grows as 



|(i?,t) = 1.2 x 10 6 e^R-^ el E dl pf t (7) 

for time t and orbital distance R. E^i is the initial dust surface density at 1 AU, and pd is 
the dust grain solid density. e s is the sticking efficiency, taken to be a constant in the range 
0.001-1. p is the radial mass exponent, where S « R~ p . The constant and exponential are 
determined semi-analytically, and have units such that the final expression is correct for cgs 
inputs. For particles above one meter, processes such as gravitational instability dominate 
over accretionary growth, and photolysis is unimportant so growth is simply turned off. 
Although eq. [7] is a simple expression and considers only a single characteristic size, it 
broadly agrees with more detailed calculations of grain growth in the micron-to-meter range 
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e.g. 



Weidenschilling 



1993). 



3.3. Photo-evaporation and Viscous Evolution 



The gas disk evolves using the standard prescription of iPringld (1l98ll ). with a = 0.01. 
For simplicity only gas is transported, not dust. Photo-evaporation by an external O 



star of brightness 



Matsuyama et al. 



5 L & is also included; the photo-evaporation model is as laid out in 



(120031 ). The photo-evaporation model considers both EUV and FUV flux 
from the external star as it removes the gas disk in 10 5 - 10 6 years. Flux from the disk's 
own central star is not included. The nominal stellar distance is 0.1 pc (giving a flux 10 6 
Go); additional cases consider distances of 0.01 pc and 0.5 pc. Photo-evaporation removes 



gas but not dust, because e xcept at the o uter edge, grains grow too quickly to be entrained 



in the escaping flow of gas ( IThroo]: 



2000). In the current model, photo-evaporation sets the 



timescale for removal of the gas disk, but does not itself affect photolysis. 

Dense clusters such as Orion also are rich with outflows, stellar winds, and close stellar 



encounters. The se effects, while easily visible in the nebu 



individual disks ( iThroop and Bally 



2008 



Allen et al. 



a, have only minimal effect on 



20071 ) . so they can be safely ignored. 



4. Results 

Simulations are presented here for our nominal well as six additional test 

cases, called Massive, High Flux, Low Flux, Slow Grow, Delay, and Debris. The initial 
conditions for each are listed in Tables [1] and |2j 

Results for the Nominal case running for 3 x 10 5 years are shown in Figs. HHSl Grains 
grow steadily and reach meter sizes at the inner edge and almost 1 cm at the outer edge. As 
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a result the opacity drops (Fig. [2]), reaching unity at the inner edge in about 10,000 years 
and dropping below this throughout the entire disk by 10 5 yr. The surface density at the 
inner edge is highest, but the rapid grain growth more than compensates for this, clearing 
this region earliest. At first, the production of organics proceeds at a uniform rate across 
the disk, because all regions of the disk are optically thick and therefore intercepting every 
single photon (Fig. [3]). As the inner-disk opacity drops, photons begin to pass through 
the disk without interacting. This causes organic production to stop in this region, and 
this 'production front' slowly moves outward in the disk until organic production has 
ceased throughout, at about 2 x 10 5 yr. Further integration beyond this point produces 
no more changes, because the disk's opacity has dropped below unity and little UV flux is 
intercepted. Looking at the total UV exposure (Fig. @|, the greatest dose is delivered at the 
outer edge, which receives ~ 3000 photons molecule -1 . The inner region receives a much 
smaller dose of <100 photons molecule -1 inward of 5 AU. The outer edge's dose is greater 
because of both slower grain growth and lower surface density. The dose of 3000 photons 
molecule -1 is be sufficient to photolyze raw ices into complex organics about seven times 
over. As seen in FigHJ the peak organic density is at ~ 20 AU. Outward of this the disk's 
raw material inventory drops, and inward the grain growth is too fast to allow for much 
photolysis. By the end of the run, total organic production is 4.5 x 10 29 g, or close to half 
of the original ice mass of 1.1 x 10 30 g (Fig. EJ). Photo-evaporation has removed the entire 
gas disk in 2.5 x 10 5 yr. Numerical results from this run and all others are in Table [2j 

Next shown (Figs. EHSJ) are results from the Delay case. In this, the disk is the same, 
but the UV source comes on only after a delay of 7 x 10 5 yr, allowing for some grain growth 
to happen before photolysis. As can be seen, the resulting organic surface density and UV 
exposure are about 10 x lower than that of the Nominal case. 

Figures [9TTTT1 show the Debris case, where the grains have been started at a constant 
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size of r = 1 m. In this case the disk's initial opacity is so low (< 0.01 at the outer edge) 
that very little UV is intercepted, and the organic production is about 10 6 x lower than 
the Nominal case. Both the Delay and Debris cases simulate different conditions for older 
disks: that is, disks which have evolved outside the HII region and then enter into it, or 
ones born before the O/B stars turn on. Low-mass stars are generally believed to form over 
the course of several Myr, before the first high-mass stars form, so most disks are probably 
somewhere between these two and the Nominal case. 

In the Slow Grow case (Figs. [E2HI3]), the grain sticking efficiency has been reduced to 
e s = 0.0001, to promote slower grain growth and longer-lived disks. As a result the total UV 
dose is increased by an order of magnitude, allowing ices to be entirely photolyzed outward 
of 7 AU. The total organic production is a few times greater than the Nominal case, and 
production would continue to increase given more time (in contrast to the Nominal model). 

Results for three additional cases are shown without figures, which are qualitatively 
similar to those presented. In the Massive run, the disk mass was increased by 10 x to 
0.4 Mq. This decreases the net UV dose per molecule in the inner disk, but in the outer 
disk where there is already more than sufficient flux in the Nominal case, the result is a 
net greater production of organics. In the High Flux case, the UV flux is increased by 
10 x to 10 7 Go, which is the highest that most disks typically encounter for brief periods. 
The effect is to increase photolysis, but by less than 10 x because the outer disk is already 
flux-saturated. Finally, in the Low Flux case, corresponding to a disk which is 0.5 pc 
distant from the cluster core, the photolysis rate is scaled downward, roughly in proportion 
to the reduced flux. 

Most of the runs share several similarities, including: a) the total raw exposure is 
sufficient to photolyze the entire initial ice abundance of much of the disk; b) the photolysis 
is greatest at the outer edge; and c) the dose in a typical cluster is high enough that in 
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10 5 yr the entire disk outward of 20 AU is photolyzed. 

The model here is simplistic in several regards. First, the laboratory yields used here 
are for photolysis of virgin ice surfaces exposed for short times, and the efficiency e p will 
drop with time as the virgin surface becomes irradiated. More processing will also result in 



the d estruction of some amino acids into nitriles and gas-phase molecules (IBernstein et al. 



2004 ). so the quantities presented here are upper limits. Heating from inward radial 
transport and shocks may cause additional destruction. The total UV dose in the Nominal 
case is high but not extreme: the highest-dose region at the outer edge receives only 7x the 
flux required for complete photolysis, and regions inside of the edge receive far less than 
this (< 1 inward of 10 AU). While eventual destruction of many of these organics must be 
assumed, the actual amount depends on parameters outside the bounds of this study. 

It is assumed that each individual grain is uniform throughout, while in reality flux is 
deposited only to the outer skin of a grain, where it can build up a crust layer. UV may 
be able to penet rate up to a few mm in r ough or fractured surfaces, but anything beyond 



this is doubtful (jOrzechowska et al 



20071 claims 1 meter penetration of UV into Europa, 
but this seems highly unlikely in the case of any realistic ice). However, most of the flux 
in the simulations is deposited when the grains are less than a few mm, making this a 
reasonable approximation. For instance, Fig. [3] shows that UV deposition has stopped by 
about 20,000 years in the Nominal disk. Fig. [1] shows that by this time, grains in this region 
remain mm-sized. In addition, grains of this size are likely to be fluffy, cracked, collisional 
aggregates which expose new surface regularly. And in all cases, the rapid grain growth 
occurs in the inner disk, where the models predict very little photolysis in the first place. 
The handling of UV deposition here is thus probably accurate enough given the model's 
other uncertainties. 

The model here assumes the disk to be made of ice only, and not silicates. In reality 
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the ice:silicate ratio is approximately 1:1 outward of the snow line, causing the net flux to 
be roughly halved for homogeneously mixed grains. Even within individual ice grains, the 
formation process may concentrate some reactants towa rd the center and others toward the 



outside, slowing photolysis further (jCollings et al. 



20031). 



In addition to chemically processing grain surfaces, UV photons can erode grains 
by photo-sputtering molecules into the gas phase, where they are easily remo ved by 



the photo-evaporative wi nd. Sputtering may be as high as 1 meter per Myr ( IThroop 



2000 



Westley et al. 



19951 ) for fresh surfaces at 10 6 Gq. However, for both photolysis and 



sputtering, large bodies will probably build up a crust of silicates which will prevent further 
effects of exposure. It is difficult to estimate this effect without more detailed laboratory 
modeling of heterogeneous mixtures of thick, collisionally processed ices and silicates. The 
interaction between sputtering and photolysis should be considered in future models. 

Our model handles the radial tran sport of gas but not d ust. Dust transport is a poorly 



Brownlee et al. 



20061 ) . and should be included in the 



understood yet important process {e.g. 
next generation of models. 

The model here finds that the disks become optically thin by 10 5 -10 6 yr, somewhat 
shorter than the observed disk lifetimes of 1-5 Myr. The difference between these two 
timescales is that our model considers just monotonic grain growth, while real disks continue 
to release small secondary dust grains during particle collisions. These continually produced 
small grains remain visible even as most of the mass of the disk grows into planetesimals. 
UV photolysis (and photo-evaporation) can therefore continue for much longer periods than 
the lower-limit timescales predicted here. The assumption of short disk lifetimes here causes 
our flux numbers to be conservative. 
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5. Discussion 

The model here leads to one robust result: there was sufficient external UV flux in the 
young solar system to allow for ice photochemistry to play an important role in the solar 
system's chemical makeup. Complex organics can be rapidly produced in flux environments 
typical of dense star clusters, especially in the outer solar system. Previous models for 
photochemistry in our solar system have not considered ice photochemistry from external 
UV sources. The current model describes a new pathway for production of complex organics 
in the young solar system. This adds to the two existing pathways of aqueous alteration 
on meteoritic parent bodies and photo-synthesis in the ISM, thus relaxing the conditions 
required in the early solar nebula. The total external UV energy available can exceed that 
of all other energy sources combined. 

Of the simplifications made here, two are worth discussing. First, photolysis is 
assumed to be a one-way process, when in reality continued UV can break apart amino 
acids that it creates. Continued laboratory work and improved modeling will constrain 
the net production rate, allowing improvements from the upper limit yields presented 
here. Second, the interaction between photo-evaporation and photolysis is simplified. 
Photo-evaporation is driven by UV absorption into dust or ice grains, and if they are 
tiny enou gh the photo-evap orative outflow can entrain these grains, as seen in some large 



outflows (jBalog et al. 



20081). Small grains might be continually produced as collisional 
debris in young disks. So, if the disks are dominated by a continually refreshed population 
of small dust grains, these grains might act as a shield, protecting the inner disk from UV 
but themselves providing a large source of irradiated organics delivered back to the ISM in 
the outflow. 

The largest free parameter in these models relates to timing: if disks evolve in a dark 
environment for more than 10 5 -10 6 yr, grains grow to meter sizes and larger where UV has 
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no chance of penetrating to cause any more than trivial amounts of photolysis. In dense 
star clusters, the formation of low-mass stars is thought to continue steadily up until the 
moment that massive O/B stars form. Once they have turned on and the molecular cloud 
is ionized, star formation stops. Low- mass stars that form just prior to the O/B stars will 
thus end up with the largest amount of organics, while the oldest disks (up to several Myr 
in the case of Orion) will have much less. 

The results here predict that organic molecules will be abundant in young circumstellar 
disks, but to date only a few of these molecules have been actually detected. However, 
help is on the way in the form of the Stratospheric Observatory For Infrared Astronomy 
(SOFIA) and the Atacama Large Millimeter Array (ALMA). SOFIA, operating in the 
3-300 /im range at a spectra resolution of ~ 100,000, will be able to detect vibrational 
and bending lines of water ices, methanol, and some CO bands. Organics such as amino 
acids are too heavy to have easily visible bands, but the structure of the water ice matrix 
as modified by these heavier molecules will be able to be studied by SOFIA. ALMA will 
provide brand new 0.1" high resolution imaging of the disks, allowing insight into the radial 
distribution of molecules across the disk. 



6. Conclusions 

• There was sufficient flux from nearby O/B stars in the Sun's birth environment to 
convert the the majority of the ice mass in the young Solar System and other proto-planetary 
disks into complex organic molecules within a few 10 5 yr. Even when considering the disks' 
high optical depth and rapid grain growth, both of which sequester surface area, enough ice 
remained exposed for photolysis to continue rapidly. 

• UV photolysis from external sources is likely to be the source of the meteoritic amino 
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acids, which have D/H ratios difficult to reconcile with the warm, aqueous environments 
required from Strecker synthesis of these species. 

• Organic production is predicted to be substantially less for disks that have evolved 
for more than 10 5 -10 6 years before O/B stars have turned on, because the grain growth in 
this time locks up surface gainst photolysis. 

• Externally driven photochemistry is most important in the outer disk, where grain 
growth is slow. The inner disk is warmer, more massive, and sees more rapid growth, and 
thus processes such as Fischer- Tropsch reactions may dominate the production of organics 
here. 

• External UV flux in 1 Myr provides roughly 100 x as much flux as from the central 
star, and 10 4 x as much energy as gravitational collapse. This energy is easily usable for 
chemical reactions because it is directly incident on ice grain surfaces. 

• Photo-evaporation can occur at the same time as photolysis of ices, but these two 
processes operate largely independently. Photo-evaporation removes only the smallest 
grains, and photolysis can continue long after the gas disk has been removed. 

Dense star clusters have historically been thought of as a hazard to planet formation, 
because the rapid timescales of disk destruction limit the conditions under which planets can 
form. Our earlier results showed that under some conditions p hoto-evaporation may in fact 



speed the formation of planetesimals (IThroop and Bally 



20051 ) . Our new results here add 



to this irony, showing that planetary systems that can form in such 'hostile' environments 
may also be be among the richest in organic molecules and and the pre-cursors of life. 
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A. Appendix: Grain Growth 



Throopl (120001 ) and iThroop et al\ (120011 ) developed a simple model for accretionary 
grain growth. The details of that model were presented only schematically, so a more 
detailed derivation is given here. In this model, particle grow by pure accretion, which is 
often considered to be the dominant process for the earliest phase of growth, where particles 
range from microns to ~ 1 meter. 

Consider a disk which has collapsed from the ISM, with an initially uniform particle 
size of ro- The disk is arranged such that its full height H at radius R is given by 



H = R/10 (Al) 

Within this disk, dust grains collide and begin to grow by simple accretion. For simplicity, 
we assume that the particles can be defined by a single grain size r(R) which is a function 
of orbital distance. A particle-in-a-box collision model gives the grain growth rate as 



dm/dt = anvme s (A2) 

where the collision cross-section is o = 2rrr 2 , n is the particle number density, v is the 
collision velocity, the particle mass is m = |7rr 3 pd, and e s = [0, 1] is the sticking efficiency 
for a single collision. Assuming vertical homogeneity and no radial transport, the particle 
volume density n can be computed directly from the disk structure as 
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n = Yi/(Hm) 

and the surface density of the dust component S d (normalized to at 1 AU), 

R 



(A3) 



mi 



IAU 



(A4) 



To compute the collision velocity v. 



Throopl (120001 ) examined the interparticle velocities 



present in the convecting, turbulent disk. The collision velocity depends on particle size 



because partic 
to the work of 



es are trapped in eddie s of different s t rengt hs. Performing semi- analytic fits 



Voelk et al. 



flf980h and 



Mizuno et al. 



(119881 ) . they approximated the collision 



velocity for dust grains in a gas disk surrounding a 1 M Q star as 



10V/ 2 (2.2 10- 16 Rf/ 4 



(A5) 



where all units are cgs. Combining these eqs. IA1HA51 integrating, and solving for r yields 
the result 



dr/dt = 1.2 x 10 6 e 42 - 6p R- 2 ~% e 2 s £ 2 p d 5/ * t 



_2 V 2 -5/3. 



(A6) 



which can then be easily used to compute r(t) or m(t). The constant and exponential 
out front are determined semi-empirically and have the proper units such that the result 
is correct for cgs inputs. The expression shows a strong exponential dependence on the 
surface density exponent p, but this is deceptive: p usually has only a small range of values 
(1-2), and the dependence is reduced by other terms so that p only has an appropriately 
small effect on the grain growth. 



Throop 



(|2000[ ) tested this equation (Fig. [T5|) and found it to agree well with accretionary 



growth calculations done numerically with full size distributions, using both their own 
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Weidenschillingl (119971 ) . In the case of 



calculations and those of iMizuno et all (119881 ) and 
e s = 1, it gives an upper limit to the grain size in a young disk growing by coagulation. For a 
standard 0.05 M & disk with gas:dust ratio of 100:1 and 0.1 /xm initial grains, eq. IA6I predicts 



growth to 
verified by 



cm in 10 4 year s and 1 m in 10 5 years inward of 30 AU. This rapid growth is 
WeidenschillingJ ( 119971 ). which finds growth of up to several meters within 10 5 



years at 30 AU. Because the model considers only accretionary growth, it is useful only 
for particles up to a few meters in size, after which gravity and settling begin to become 



important. It also ignores radia 



such as gravitational instability (IJohansen et al 



drift, fragmentation , and possibly faster mechanisms 



20071 ). Nevertheless, for computing rough 



particle size estimates and their dependence on disk parameters in the youngest disks, it 
provides useful limits which are consistent with more complex codes. 
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Parameter 


Value, Nominal 


Surface Density 


v r->— 3/9 

£ oc R 6/2 


Initial grain size 


r = 0.2 fj,m 


Photolysis efficiency 


e p = 0.0025 


Sticking efficiency 


e s = 1.0 


Stellar mass 


M = 1 M . 


Initial disk mass 


M d = 0.04 M 


Disk edges 


T~t f\ A A T T 7 1 o r\ ATT 

Ro = 0.4 AU; Ri = 30 AU 


Ice: Gas mass ratio 


1:100 


UV flux 


F = 10 6 G 


UV start time 


^uv = yr 


Run time 


W = 3 10 5 yr 



Table 1: List of parameters, Nominal case. 
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Trial Name 


Changed Parameter 


Organic Production 


Nomi 1 

l\l Will _L J.-L CJ- _L 




4.5 10 29 g 


Massive 


M od 4O.4M 


1.6 10 29 g 


High Flux 


F -> 10 7 G 


7.0 10 29 g 


Low Flux 


F -> 2.5 x 10 4 G 


3.1 10 28 g 


Slow Grow 


e s -> 0.001 


6.3 10 29 g 


Delay 


tuv -> 700,000 yr 


2.9 10 28 g 


Debris 


ro — )■ 1 m 


4.0 10 23 g 



Table 2: List of parameters varied. 



Source 


Flux at 10 AU 


Reference 


UV flux inside dense molecular cloud 


10 -4_ 10 -5 Gq 


Prasad and Tarafdar (1983) 


Present day, interstellar 


1 Go 


Habiner (1968) 


Present day, Solar 


10 G 


Cox (2000) 


Young Sun, reflected from IPM 


1000 G 


Gladstone (1993) 


Central T Tauri onto flared disk 


3000 G 


Alexander and Armitaee (2007) 


External stars, small cluster 


10 3 10 5 G 


Adams et al. (2006) 


External stars, large cluster 


10 5 -10 7 Go 


Johnstone et al. (1998) 



Table 3: UV fluxes in different environments. 
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Fig. 1. — Grain size, Nominal. Grains grow to almost cm sizes throughout the disk within 
3 x 10 5 yr. Growth is stopped at 5 m (black) as other growth mechanisms begin to take 
over. 




Fig. 2. — Opacity, Nominal . By the end, the opacity has dropped to below 1 throughout 
the disk as grains grow. In the regions of low opacity, photolysis is inhibited because the 
disk intercepts little flux. 
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Fig. 3. — Organic production, Nominal. The integrated flux and opacity is converted into a 
maximum organic photolysis yield. The peak density is reached near the outer edge: inward 
of this grains grow rapidly, and outward the disk density drops off. 
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Fig. 4. — Photons per molecule, Nominal . This shows the net exposure of each original ice 
molecule in the disk to UV flux. Grains at the outer edge receive the highest flux because 
of their slow growth and the low surface density. 
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Fig. 5. — Mass evolution, Nominal . As ices are converted into organics, the solid and dashed 
lines approach each other. The dotted line plots the mass of the gas disk, as it is dispersed 
by photo-evaporation and viscous evolution. 
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Fig. 6. — Photolysis, Delay. In this case the disk is given a head-start of 700,000 years to 
evolve before the UV flux is started. Grains grow, and the total organic production is about 
5% of that in the Nominal case. 




Fig. 7. — Photons per molecule, Delay . 
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Fig. 8. — Mass evolution, Delay. 
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Fig. 9. — Photolysis, Debris . Here the grains start as a uniform 1 m size, with the same 
radial distribution as in the Nominal case. Organic exposure is very low, about 10 -6 that of 
the Nominal model. 




Fig. 10. — Photons per molecule, Debris . 
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Fig. 11. — Mass evolution, Debris. 
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Fig. 12. — Photolysis, Slow Grow. The grain sticking coefficient is reduced from 1.0 to 10~ 3 . 
Grain growth is slowed, and exposure to UV increases as a result because of long-lasting small 
grains. 




Fig. 13. — Photons per molecule, Slow Grow. 
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Fig. 14. — Mass evolution, Slow Grow. In this case more than 50% of the original ices can 
be photolyzed into organic molecules. 
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Fig. 15. — Grain size evolution. Plot show accretionary grain growth according to eq. IA6| 
for disk of p = 3/2 and 0.1 M ; growth is stopped above 1 meter radius. Radial bins are 
spaced logarithmically. 



